Microscopic description of dissipative dynamics of a level crossing transition 
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We analyze the effect of a dissipative bosonic environment on the Landau-Zener-Stiickelberg- 
Majorana (LZSM) level crossing model by using a microscopic approach to derive the relevant 
master equation. For an environment at zero temperature and weak dissipation our microscopic 
approach confirms the independence of the survival probability on the decay rate that has been 
predicted earlier by the simple phenomenological LZSM model. For strong decay the microscopic 
approach predicts a notable increase of the survival probability, which signals dynamical decoupling 
of the initial state. Unlike the phenomenological model our approach makes it possible to study the 
dependence of the system dynamics on the temperature of the environment. In the limit of very 
high temperature we find that the dynamics is characterized by a very strong dynamical decoupling 
of the initial state — temperature- induced quantum Zeno effect. 
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I. INTRODUCTION 

The exactly soluble Landau-Zener-Stuckelbcrg- 
Majorana (LZSM) model [l|-[4| is a very popular tool 
for estimating the transition probability between two 
quantum states whose energies cross in time. A level 
crossing in combination with adiabatic evolution is 
of great practical significance for it leads to complete 
population transfer between the two states; various 
experimental level crossing techniques have been devel- 
oped and demonstrated over the years @. Insofar as a 
quantum system is often immersed in a noisy environ- 
ment, a number of studies have dealt with the effects of 
dissipative environments on this model [(HI, [HI [Hj] . 

Despite its extreme simplicity — linearly changing en- 
ergies and a constant coupling — the LZSM model often 
provides unexpectedly accurate results when applied to 
real physical systems with elaborate time dependences. 
This feature, and the appealing simplicity of the tran- 
sition probability, have made this model vastly popular. 
There are several additional intriguing features associ- 
ated with the LZSM model. It turns out that some 
standard approximate methods, e.g. the Dykhne-Davis- 
Pechukas approximation [l3j , when applied to this model 
give the exact transition probability. Moreover, when the 
upper state is allowed to decay irreversibly to other levels 
outside the two-state system with a decay rate T the sur- 
vival probability, i.e. the population of the lower state, 
does not depend on T 0, Q , while the population of the 
upper decaying state vanishes at large times. This is sur- 
prising because the evolution of the ground-state popula- 
tion does depend on T: it approaches its asymptotic value 
smoothly for nonzero T and in an oscillatory manner for 
r = 0. It is also remarkable that an adiabatic elimina- 
tion of the decaying state provides the exact result for 
the ground-state population. It has been demonstrated 
that all these features are unique for the LZSM model 
and do not apply to any other model Q. 

The latter results for the lossy LZSM model have been 
obtained in the framework of the phenomenological in- 



clusion of the decay rate in the Schrodingcr equation 
as an imaginary energy term in the Hamiltonian 0, [1] . 
This approach is known to be justified for weak decay. 
We have demonstrated recently that a rigorous micro- 
scopic master equation treatment, developed in Ref. flij ]. 
can deliver dramatically different results compared to 
the phenomenological approach in the description of the 
stimulated Raman adiabatic passage technique in open 
three-state quantum systems in the regime of strong de- 
cay [HI, [l6| . It is therefore interesting to examine the be- 
havior of the LZSM model in an open quantum system by 
the master equation approach and verify its properties, 
e.g., the validity of the independence of the survival prob- 
ability on the decay rate. Moreover, the microscopic ap- 
proach allows us to study the LZSM model in the regime 
of nonzero bath temperature, which is inaccessible by 
the phenomenological approach. In view of the numer- 
ous applications of the LZSM model and the fact that a 
quantum system is often subjected to decoherence, our 
results are of potential significance in a number ofphys- 
ical systems 



17 1, such as J osep hson junctions [l8|, Qi 
cold atoms in optical lattices [2(J , spinorial Bose-Einstein 



condensates 1211, and others 



The paper is structured as follows. In the next sec- 
tion we present the physical system, by starting from 
the ideal case and then introducing the interaction with 
a bosonic environment and deriving the relevant mas- 
ter equation. In the third section we compare numer- 
ical simulations obtained by our microscopic approach 
with predictions coming from the previously studied phe- 
nomenological model. Furthermore, the dependence of 
the dynamics on the temperature of the environment is 
simulated and explained by invoking the notion of dy- 
namical decoupling. Finally, in the last section, we give 
some conclusive remarks. 
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II. PHYSICAL SYSTEM 



B. The microscopic dissipative model 



A. The Ideal Case 

The ideal model comprises two quantum states |1) 
|2) coupled by the time-dependent Hamiltonian 



H(t) 



A(t) 



and 



(1) 



In order to model the interaction of this system with a 
dissipative environment, we add a third state |3). There- 
fore we consider a three-state system described by the 
Hamiltonian (h = 1) 



Hs(t) 



-A/2 n 
Q A/2 
-u> 3 



(2) 



We assume this model to be "exact" : this is possible for 
instance if the two states |1) and |2) correspond to two 
Zeeman sublevels of an atom or a molecule splitted by 
a time-dependent magnetic field, while the coupling Q 
is induced by a static electric field in the presence of a 
dipolc operator which has only real matrix elements, i.e., 
d(\l) (2| + \2) (1|). It is worth noting that, in the deriva- 
tion of the Hamiltonian for such a system, no rotating 
wave approximation is performed. 

The eigenvalues of Hg are ±e, with e = y/Q 2 + A 2 /4, 
and uj 3 . The corresponding cigenstates are 



|+) = cos^ |1) + s'mtp |2) , 
|— ) = — sin <p |1) + cos (f |2) 

13), 



(3a) 
(3b) 
(3c) 



where tan cp = (A/2 + e)/0. The quantities A and fl are, 
in general, time-dependent. In the following we will focus 
on the case = (i.e., static electric field) and A = n 2 t 
(i.e., linearly changing magnetic field), with t spanning 
the range [— t, t]. 

According to the calculations made in [H-Q| , where the 
time of interaction r is assumed to be infinite, when the 
system is prepared in state |1) at t = —00, the survival 
probability of state |1) at t = 00 is given by 



Pi (00) 



(4) 



In the adiabatic limit, Q/k 3> 1, the population of the 
initial state |1) vanishes, Pi (00) — > 0, and complete pop- 
ulation transfer |1) — > |2) takes place. This formula is 
approximately valid also for finite values of r provided 
t 3> fl/n 2 [9J, because the coupling il induces transi- 
tions between |1) and |2) only in the proximity of the level 
crossing. In other words, the LZSM formula ^ applies 
to finite times t which are much larger than the char- 
acteristic transition time, which in the nearly-adiabatic 
regime is of the order of Vl/n 2 0. The values of r con- 
sidered in the numerical simulations below will always 
satisfy the above condition. 



Let us now consider the interaction with a bosonic 
bath, so that the system is described by 

H = H S + H B +H ll (5a) 

where 

H B = ^2 0Jka{a k , (5b) 

k 

£ri = (|2><3| + |3)(2|)®£ 5jk (a fc + 4) . (5c) 

k 

The interaction Hamiltonian in (I5cl) has the form 



Hi = A ® B , 

A =|2> (3| + |3) <2| 
B = ^fjk (a k + a\ 



(6a) 

with 

(6b) 
(6c) 

Following the general theory of Davies |l4| and the rele- 
vant consolidated approach [l5|, [lj| [22[ , the master equa- 
tion can be written down, by using standard methods 
[23l l24j , in terms of the instantaneous jump operators 
between the eigenstates ([3]) of the system Hamiltonian, 

(7a) 
(7b) 
(7c) 
(7d) 
(7e) 



0. 



A(w + _) = A 3 (cj_+) 
A(lj +3 ) =smip\3) (+| , 
A(u 3+ ) = sin^|+) (3| , 
A(uj_ 3 ) = cos</?|3) (-| , 
A(w 3 -) = cos^|-) (3| . 



Moreover, since we are in the presence of a level cross- 
ing we may not be allowed to perform the rotating-wave 
approximation. Therefore the master equation in the 
Schrodinger picture is 

-At(w / )A(w)p]+H.c, (8) 



with 



r(«) 



s tr B [B\t)B(t - S ) PB (Q)] , (9) 



where B(t) is the operator B in the interaction picture 
at time t. In the continuum limit we have 



r(u) = | 5 (M)| 2 D(M)(i + iVM) 



(10) 



where g(co) is the coupling constant g k corresponding to 
the mode of frequency w, D(uj) is the density of the bath 
modes at frequency oj, and N(lo) is the number of pho- 
tons in each mode at frequency u. Since the latter quan- 
tity is the occupation number for a thermalized boson 
system: N{u) = l/(e" /e - 1), where 6 = k B T with T 
being the temperature and k B the Boltzmann constant, 
the sign of ui is automatically taken into account since 
when u> < one has 1 + N(u) = -N(\lj\). 
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FIG. 1: (Color online). Population of the state |1) at t = r 
vs the decay rate V (in units fi) at O = 0, for different values 
of the chirp parameter k: k/CI = 1 (red solid line), k/Q = 2 
(green dashed line), k/Q. = 4 (blue long dashed line). Here 
we have assumed uj^/fl — 10 3 . Each time evolution occurs in 
the time interval [— r, r], with r = 30fi - . The initial state 
of the system is p(—r) = |1) (1|. 



III. ASYMPTOTIC POPULATION OF THE 
GROUND STATE 



A. General behavior 



FIG. 2: (Color online). Population of the state |1) at t = t 
as a function of the duration r (in units of f2 _1 , spanning the 
range [10f2 _1 ,60n -1 ]) and of the decay rate T (in units of 
Q). The dependence of the final population on the decay rate 
is not affected much by the duration. 




We have integrated numerically the master equation 
(|SJ) in order to find the behavior of the survival proba- 
bility of the state |1) after the interaction. We consider 
for simplicity a fiat spectrum for the reservoir, i.e. we 
assume |g(M)| 2 -D(|u;|) = T independently of u. 

In Ref. Q a phenomenological model describing the 
level crossing transition in the presence of losses toward 
external states has been studied by adding an imaginary 
part (— iT) to the energy of the state |2). The analy- 
sis there and previous work Q showed that the survival 
probability in the LZSM model does not depend on the 
decay rate at all. However, in the finite version of the 
LZSM model where the interaction duration is finite @, 
and in several other analytic models, the phenomenolog- 
ical approach predicts an increase of the survival proba- 
bility in the presence of strong damping (dynamical de- 
coupling): Pi(oc) — > 1 when T — > oo. 

Our approach gives essentially the same zero- 
temperature dynamics as the one predicted by the phe- 
nomenological model for the finite LZSM transitions. 
Figure [T] shows P\{t) (for r J> fi/ft 2 ) as a function of the 
decay rate T, for different values of the chirp parameter 
k. As k increases, the adiabaticity condition !1/k>1 de- 
teriorates and more population is left in the initial state 
at T = 0. The figure clearly demonstrates that for low 
decay rates the survival probability is independent of T, 
while it increases for larger values of T. In the strong 



FIG. 3: (Color online). Population of the state |1) at t = r vs 
9 (in units fi _1 ), for T/fi = 0.1 (green solid line), T/Q = 1 
(red dotted line), F/Q, = 25 (blue dashed line), T/Q, = 100 
(black long dashed line). Here we have k/Q = 1, ojs/Q — 10 s . 
Each time evolution occurs in the time interval [— r, r], with 
t = 30Q~ . The initial state of the system is p(—r) = |1) (1|. 

damping limit the survival probability approaches unity 
as in the phenomenological model. 

Figure [5] shows the dependence of -Pi(t) on both T and 
r, revealing that the dependence of the final population 
on the decay rate is only slightly affected by the duration 
time r, and that for smaller values of T the influence of 
t is even less visible. However, we note that, increasing 
the value of t, there is a slight shift of the region where 
the final population starts increasing with respect to T. 
Since the larger r the larger the value of T needed to see 
a significant increasing of the population Pi(t), one can 
conjecture that in the limit t — > 00 one can recover the 
independence of the population from r predicted by the 
phenomenological model. Nevertheless, developing this 
analysis is beyond our present numerical resources. 

Differently from the phenomenological approach which 
works at zero temperature, our approach is able to pre- 
dict the behavior of the system also at nonzero temper- 
ature. Figures [3] and 0] show the survival probability as 
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a function of the temperature for different values of T. 
The two figures refer to different values of the parameter 
k: k = 1 and k — 2 respectively. In both figures, one 
can see that there is a low temperature range where the 
behavior at zero temperature is reproduced. Increasing 
the temperature, there is a zone with a more complicated 
dynamics, and in all cases one can see that for very high 
temperature the survival probability approaches unity. 
This is a universal feature of the dynamics related to this 
dissipative LZSM model, as we will show in the following. 



FIG. 4: (Color online). Population of the state |1) at t = r vs 
6 (in units fi -1 ), for r/fi = 0.1 (green solid line), r/fi = 1 
(red dotted line), r/fi = 25 (blue dashed line), r/fi = 100 
(black long dashed line). Here we have «/fi = 2, uj/fl = 10 3 . 
Each time evolution occurs in the time interval [— r, r], with 
t = 30fi~ . The initial state of the system is p(— t) = |1) (1|. 



B. Explanation of the behavior at high 
Temperature 

In the high-temperature limit we have N(lo) w 0/w, 
so that the rate equations can be written down as follows: 



Pn 

P\2 
Pl3 

hi 

P23 
P33 



TQesm2ip\ 



ifl 



H^3 + y ) 



Pn + 



r0(ecos2iy9 — lu 3 ) 



TQesm2ip 

l\lp22 H % 2~ P33 : 



r9(ecos 2<p — ljs) 



2(e 2 -^ 2 ) 



Pl3 - iQp23 + 



-id 
2T ( 1 



TGe sin 2if 



Pl2 



6 



uis — e cos 2ip 



iCl 



P33 ■ 



TQe sin2ip 



P21 



P12 

r6esin2y 
— 2 T~ P32 

2re (-w 3 + ecos2< /3 ) 



= -in - 



T6e sin 2tp 



P13 



+ r (^1 + 20 

r0esin2(yj 



u>3 — e cos 2ip 



1(^3 - y) - r 



P32 , 



TQ(e cos 2(p — 0J3) 



P23 



-P22 



T6e sin 2tp 



-P31 



r9esin2(^ 2r9 (-L03 
2~P 12 + ~2 — 772~P 21 72~ 



€ cos 2(p) 



P22 - 2T 1 + e 



LJ3 — e cos 2(p 



P33 



(11a) 
(lib) 

(lie) 



(lid) 

(lie) 
(llf) 



and P21 = p*2, P31 = P13, P32 = f>23- Let us introduce the perturbation parameters r\ = W3/O <C 1 and £ = e/w3 <C 1 
(recall that f2 < e and A < e) and define the symbol O x = O(cc), i.e. infinitesimal of the order of x or smaller. On 
this basis, after introducing the vector notation for the density operator p = (pn, P12, P13, P21, p22-> P23, p3i, P32, P33), 
one can rewrite the rate equations (|lla[) - (|llf|) as dp/dt = Cp, with 
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where we have used such equivalences as fl/Q = fl/co^ x UJ3/Q — 0(£rj), etc. Keeping only first order terms in £ and 
77, one obtains: 
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from which we conclude that the projector |1) (1| is 
an eigenoperator of the Lindbladian corresponding to 
the eigenvalue 0, at first order in rj and £. This ex- 
plains why the population of the state |1) is preserved 
at high temperature. In fact, we are in the pres- 
ence of a temperature-induced dynamical decoupling (or 
temperature- induced quantum Zcno effect). Though this 
effect resembles the dynamical decoupling induced by de- 
cay (25], [2(| , it is a quite different phenomenon since it 
derives physically from the presence of thermal photons 
and the ensuing growth of decay and pumping rates at 
very high temperature. 

It is important to note that the factor G does affect 
the eigenvalues of the Lindbladian, but it does not affect 
the structure of the eigenstates. 

IV. CONCLUSIONS 

In this paper we have analyzed the Landau-Zener- 
Stiickclbcrg-Majorana model in the presence of a dis- 
sipativc bosonic enviroment by using a microscopic ap- 
proach to derive the relevant master equation. Because 
of the presence of a level crossing, the rotating- wave ap- 
proximation cannot be made for the system-environment 
coupling. The time evolution predicted by our micro- 
scopic model at zero temperature and weak dissipation 
is nearly identical to the time evolution predicted by the 
well-known phenomenological approach or the original 
LZSM model (assuming a constant coupling of infinite 



duration) Q: the survival probability of state |1) is in- 
dependent of the decay rate. For strong dissipation we 
find that the initial-state population approaches unity, as 
in the finite LZSM model [9| and other analytic models 
with couplings of finite duration . 

An important advantage of our microscopic approach 
is that it can describe also the dissipative dynamics for 
an environment at a nonzero temperature. A distinct 
feature in this case is the survival probability approach 
to unity for very high temperatures. The physical origin 
of this effect is the dynamical decoupling between state 
|1) and the subspace spanned by states |2) and |3). This 
is induced by the very strong coupling between states |2) 
and 1 3) mediated by the interaction with the environ- 
ment, which becomes stronger as the number of photons 
increases. It is important to note that in order to verify 
these results the initial state (|1)) must not be electro- 
magnetically coupled to the external state (|3V), so that 
all dissipation is due to state |2), as in Ref. 0. There- 
fore this very peculiar feature of the high-temperature 
dynamics should be measurable in experimental situa- 
tions where only one of the two states is unstable. 
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